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Abstract: When a quantum system is divided into subsystems, their entanglement entropies 
are subject to an inequality known as strong subadditivity. For a field theory this inequality 
can be stated as follows: given any two regions of space A and B, S(A) + S(B) > S(AUB) + 
S(A D B). Recently, a method has been found for computing entanglement entropies in any 
field theory for which there is a holographically dual gravity theory. In this note we give a 
simple geometrical proof of strong subadditivity employing this holographic prescription. 
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Entanglement entropy is an important tool in the study of quantum information (see 
e.g. Q). It quantifies the extent to which the state of a given subsystem of a quantum 
system is correlated with that of the rest of the system. Entanglement entropy enjoys a 
crucial mathematical property called strong subadditivity |2], ^] . Recently, Ryu and one of the 
authors of the present paper [||, ||] proposed a relationship, applicable to any quantum field 
theory with a holographic gravity dual, between the entanglement entropy of a region of space 
in the QFT and the area of a certain minimal surface in the dual spacetime. An important 
test of the validity of this proposal is whether it has the property of strong subadditivity. This 
question was investigated in a variety of examples in the paper Q , always with an affirmative 
answer. In this note we give a general argument that it does, based only on general properties 
of holographic dualities. Besides giving support to the proposal, our argument provides 
an intuitive, geometrical way to understand strong subadditivity, a property whose formal 
algebraic proof is highly non-trivial. 

The same argument can be used to establish a concavity property for holographically- 
computed Wilson loop expectation values when the loops are coplanar. This is briefly dis- 
cussed in the last section of the paper. 

1. Strong subadditivity 

The von Neumann entropy of a density matrix p, 

S(p) = -Tr(plnp), (1.1) 

quantifies the extent to which the state represented by p fails to be a pure state. If p is 
obtained by tracing over part of the Hilbert space representing a subsystem — for example, 
one that is inaccessible to the experimentalist — then S(p) is referred to as the entanglement 
entropy of the remaining subsystem. More formally, if the Hilbert space of the full system 
factorizes into Hilbert spaces of two subsystems, TC{ u ii = T~ii <8> Ti-2, then for each subsystem 
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we define a reduced density matrix p\ = Tr^ 2 pf u n, p 2 = Tr^ pf n \\, and a corresponding 
entanglement entropy S(p\) and S(p 2 ). On the basis of the concavity of the function — xlnx 
and elementary properties of Hilbert spaces, these can be shown quite generally to obey the 
following inequalities, 

\S(pi) - S( P2 )\ < S(pmi) < 5(pi) + S(p 2 ) . (1.2) 

This property of entanglement entropy is known as subadditivity (the first inequality is also 
called the Araki-Lieb inequality j7j). In particular, if the full system is in a pure state then 
the two subsystems have the same entanglement entropy. 

Now suppose the system is made up of more than two subsystems, TL{ U \\ = Hi- Then 
the inequalities (1.2) can be strengthened to yield (2|, [j| 



S{pi2) + S( P23 ) > S( P2 ) + S(p 123 ) , S(p 12 ) + S(p 23 )>S( Pl ) + S(p 3 ), (1.3) 

where p± 2 is the reduced density matrix for TCi®Tt 2 , etc. These two inequalities can be shown 
to be equivalent by the formal device of adding a fourth subsytem such that /O1234 is a pure 
state. This property of entanglement entropy is known as strong subadditivity, and its proof 
is highly non-trivial (although again it depends only on elementary properties of Hilbert 
spaces). 1 Strong subadditivity represents the concavity of the von Neumann entropy and 
is a sufficiently strong property that it essentially uniquely characterizes the von Neumann 
entropy |]| [b§ [D|. 



In the context of a quantum field theory, a natural type of subsystem to consider is that 
associated with a given region of space. To any region A is associated a Hilbert space Ha, and 
for two disjoint regions A and B we have TLavjb = ~Ha ®1~Lb- For the entanglement entropy 
associated to TCa we will write simply S(A) (rather than S(pa))- Due to the infinite number 
of degrees of freedom involved in a field theory, S(A) typically suffers from an ultraviolet 
divergence proportional to the surface area of A ]16| , |17| . In order to deal with finite quantities 
one must impose a UV cutoff (and, if the surface area of A is infinite, an IR cutoff as well). 
One may also consider subtracted quantities that remain finite as the UV cutoff is removed, 
such as the mutual information, 

I(A,B) = S(A) + S(B)-S(AUB), (1.4) 

defined when A and B (and their surfaces) are disjoint. By (1.2) this is non-negative. By 



employing these quantities, Casini and Huerta [18] showed that an analogue of the c-theorem 
in two dimensional QFTs can be derived from strong subadditivity. 

To avoid confusion, it is important to remember that the concept of entanglement entropy 
refers to a specific state of the system at a specific time. Therefore all of the regions and 
surfaces we consider in this paper are restricted to a fixed constant-time slice of the field 
theory's spacetime. 
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Figure 1: A constant-time slice of a spacetime on which a gravity theory lives, and the conformal 
boundary on which its holographically dual field theory lives. A is a region of the boundary; m is the 
minimal hypersurface in the bulk ending on dA; and r is a region of the bulk such that dr — A U m. 



2. Holographic entanglement entropy 

Recently, a proposal has been made in ^ for how to compute the entanglement entropy 
of a region of space in any quantum field theory that admits a holographic gravity dual. 
The proposal is very simple. The gravity theory lives in a space which as usual we call the 
bulk, and the QFT on its conformal boundary. (To avoid confusion, we will reserve the 
term "boundary" for the space on which the QFT lives, and use the term "surface" for the 
boundaries of various regions in the bulk and boundary.) We consider all hyper surf aces 2 m 
in the bulk that end on dA, and ask for the one with minimal area. (See Figure |].) We then 
have 

4LrN rn:am=oA 

where a(m) is the area of m. For the case when the bulk gravity theory lives on a static 



asymptotically AdS spacetime, Fursaev |19[ has given a derivation of (|2.1|) using Euclidean 



Quantum Gravity and the basic principles [£C| £lj of the AdS/CFT correspondence |22|]. 
Notice that the expression ( p.l[ ) coincides with the Bekenstein-Hawking formula of black 
hole entropy if we replace the minimal surface with a black hole horizon. Indeed, at high 
temperature the spacetime of the gravity theory generally includes a horizon; when part of 



Alternative proofs, pedagogical expositions, and reviews can be found in [j], [| [uj, |ll|, Q 
2 We use 
dimension 2 



2 We use the term "hypersurface" because m is spatially co-dimension 1; in spacetime m is of course co- 
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the minimal surface wraps the horizon, its contribution corresponds to the usual thermal 
entropy. We will see an example of this situation when we come to Figure || below. 

Three refinements should be made to (2.1). First, both sides are divergent; the left-hand 
side is ultraviolet divergent as discussed above, while the right-hand side is infrared divergent 
due to the infinite proper distance from any point in the bulk to the conformal boundary. 
It is easy to see that the latter divergence, like the former, is proportional to the surface 
area of A. In fact, these two divergences are the same, a manifestation of the usual UV/IR 
correspondence characteristic of holographic dualities. Therefore ( |2.1| ) is meant to apply in 
the presence of a UV cutoff in the QFT, corresponding to an IR cutoff in the gravity theory. 
The simplest such cutoff is a brute force one that cuts off the bulk space at a finite value of 
the holographic coordinate. The exact choice of cutoff will not be important in what we say 
below, and for simplicity of presentation we will leave it implicit in the discussion. 




Figure 2: Two examples in which the QFT lives on a compact space and the bulk contains a black 
hole. (Technically, since we are considering an eternal black hole in static coordinates, in each case the 
full spacetime consists of two copies of the region shown connected by an Einstein-Rosen bridge; this 
won't affect our discussion.) The boundary is divided into two regions A and B. Since dA = dB, if the 
bulk had trivial (d — l)st homology the corresponding minimal hypersurfaces tua and mj would be 
identical, and we would have S(A) = S(B). However, due to the requirement that each hypersurface 
be homologous to the corresponding boundary region, raj can either (left) wrap around the other 
side of the event horizon, or (right) separate into two connected components, one being tua and the 
other the event horizon. In the latter case we have S(B) = S(A) + 5bh, where S*bh is the black hole's 
Bekenstein-Hawking entropy; since Sf u n — Sbh, the Araki-Lieb inequality is saturated. 

Second, there is a complication that occurs when the bulk has non-trivial {d — l)st 
homology (where d is the spatial dimension of the bulk, which is also the spacetime dimension 
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of the boundary). This will be the case, for example, when the bulk contains a black hole. 



Fursaev's derivation of (2.1) then tells us that we should minimize a not over all hypersurfaces 
ending on dA but only over those that are homologous to A; that is, there should exist a 
region r of the bulk such that dr = AUm. See Figures [l] and ^ for examples. (See [Q, || IS] 



for further discussion.) This rule will be essential in what follows. 

Third, the formula ( |2.1| ) is exact in the limit that the gravity in the bulk is controlled 
by the Einstein-Hilbert action. Higher curvature corrections to the bulk action will lead to 
corrections to the functional a(m). For example, Fursaev |19| ] showed that, if the bulk action 
is corrected by a Gauss-Bonnet term, then a(m) is corrected by an Einstein-Hilbert term, 



a{m) 



I s/h(l + 2aR{h))+Aa I , (2.2) 

J m J dm 



where h is the induced metric on m and a is the coefficient of the Gauss-Bonnet term in the 
bulk action (see flSj for details). In order to make the variational problem for m well defined, 
we have also included a Gibbons-Hawking boundary term; 7 is the induced metric on dm, 
and K is the trace of its extrinsic curvature (in m). 

All of these regions and surfaces — both on the boundary and in the bulk — must lie on 
a single constant-time slice. In order to have a well-defined notion of "constant-time slice" 
in the bulk, we must restrict ourselves to states for which the bulk geometry is static. A 
covariant generalization of ( [2.1| ) to time-dependent geometries will be discussed in p3fl . We 
leave the proof of strong subadditivity in that context to future work. 

3. Proof 



In the paper Q the authors investigated in a variety of examples whether the formula (2.1) 



for the entanglement entropy satisfied the property of strong subadditivity, and in all cases 
studied it did. Here we will give a simple argument that it does in general. 
We begin by rewriting the inequalities fll.3j) in the forms 



S{A) + S(B) > S(A U B) + S(A n B) , S(A) + S(B) > S(A\ B) + S(B \ A) , (3.1) 

where A \ B = A n B c . We will prove the first inequality; the proof of the second one is very 
similar and is left as an exercise to the reader. 

Let niA, ms be the minimal hypersurfaces in the bulk ending on dA, dB respectively, 
and rA, fs be the corresponding regions of the bulk (so that dr^ = AUrriA, drs = BUvtib)- 
(See Figure ||, left side.) We now define the regions 

r AuB =r A Ur B , r A nB = r A r\r B . (3.2) 

We can decompose the surfaces of these regions as usual into a part on the boundary and a 
part in the bulk, 

dr A uB = {A U B) U m A uB, dr A nB = (A n B) U m An B ■ (3.3) 
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Figure 3: Two overlapping regions A and B of the boundary, with (left) their respective minimal 
bulk hypersurfaces m,A, ms and bulk regions ra, re, and (right) their minimal hypersurfaces uia and 
m,B cut up and rearranged into two new hypersurfaces juaub (the bulk part of the surface of ta U tb) 
and rriAnB (the bulk part of the surface of ta H rg)< m^yg and rriAnB end on 9(AUB) and d(A n B) 
respectively (although they are not necessarily the minimal such hypersurfaces). 

(See Figure ||[ right side.) Clearly tuaub ends on d(A U B). While nothing says that it is 
the minimal hypersurface ending on d(A U B), its area is an upper bound on the area of the 
minimal one, and therefore on 4Gn S(A U B)\ similarly for An B. Now the hypersurfaces 
n^AuB and ttiadb ar e simply rearrangements of itia and mg (meaning that ttiaub U rriAnB = 
U mfl), so they have the same total areas, 3 

a(rn A uB) + a(m A nB) = a(m A ) + a(m B ) , (3.4) 

which completes the proof. 

Note that equation ( |3.4| ) holds not just if a is the area, but if it is any extensive functional 
of the hypersurface. This means that if m and m! are two disjoint hypersurfaces with a 
common boundary, dm Pi dm' ^ 0, then we have a(m U m') = a(m) + a{m'). This is true, for 
example, for the Einstein-Hilbert term (with boundary term) added in ( ^.2[) . 

4. Discussion 

In this letter we gave a simple geometric proof of strong subadditivity of entanglement entropy 
based on the holographic formula (|2.1| ) . The extra dimension in the holographic dual obviously 
plays an essential role in this proof. Since the strong subadditivity of entanglement entropy 
should be true in any quantum mechanical many-body system, our result shows that the idea 
of holography is consistent with any quantum system from this basic viewpoint. 

3 In this sentence we've assumed the generic situation that vtia and ms intersect along (spatially) codi- 
mension 2 submanifolds. More generally we have tuaub U mAns C tua U tub and a(m.4us) + a(?riAns) < 
a{rriA) + a(ms). 



It is interesting to ask when the inequalities ( |1.3| ) are saturated. The only examples we 
know in the holographic context involve only two disjoint regions, and therefore reduce to the 
saturation of weak subadditivity, inequalities (|F^) . (It would be interesting to find examples 
where this is not the case.) The first of these, the Araki-Lieb inequality, is obviously saturated 
when the full system is pure; then each entanglement entropy is due only to correlations 
between the subsystems, rather than to the full system being in a mixed state. A system 
that is in a mixed state but nonetheless saturates the Araki-Lieb inequality is depicted on the 
right side of Figure The fact that mg is disconnected suggests that here the entanglement 
entropy of B has two separate and unrelated origins: the thermal entropy of the full system 
(jSfuii), and the correlations between A and B (S(A)). On the left side of that figure, where 
mg is connected, the inequality is not saturated. 

As for the second inequality in ( |1.2| ), it is saturated (i.e. the mutual information vanishes) 
when two regions are sufficiently far apart that their union's minimal hypersurface does not 
connect them but instead is simply the union of their respective minimal hyper surf aces. This 
was seen in explicit examples in Q. The mutual information vanishes if and only if the two 
systems are uncorrelated, i.e. p\2 = pi ® pi [fjj. It is interesting that the correlations can go 
strictly to zero in a field theory (in the large N limit). 

Finally, it is useful to notice that our argument can be directly applied to the holographic 
derivation of a concavity property of coplanar Wilson loops ||, which is closely related to 



the Bachas inequality [25]. If the curves C A = dA and C B = dB lie in the same two- 
dimensional plane, then it is clear that the holographically computed expectation values of 
the corresponding Wilson loops satisfy 

(W(C A ))(W(C B )) < (W(C AnB ))(W(C AuB )) , (4.1) 
(W(C A ))(W(C B )) < (W(C A \ B ))(W(C B \ A )), (4.2) 



where we defined C A nB = d{A n B) etc. They are equivalent to fl3.1| ) once we remember that 
the holographic Wilson loop expectations can also be found from the minimal surface [26, 27]. 
The evidence from the gauge theory side for these relations will be discussed in [28]. 
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